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Abstract. Using a deprojection technique, we study the X-ray properties of the galaxy cluster Abell 1835 observed 
with Chandra, including temperature, abundance, electron density, gas mass fraction, and total mass. A compar- 
ison with the results without deprojection shows that the properties do not change much. When we compare the 
results with those of XMM-Newton, the difference between the temperature profiles derived from Chandra and 
XMM-Newton data still exists, even if the point-spread function effect of XMM-Newton is corrected. To inves- 
tigate the reasons for the difference, we used the spectra to cross-calibrate the temperatures. They show that 
the Chandra spectra can be fitted well with XMM-Newton temperatures. Furthermore, we derive the electron 
density profile from Chandra data with XMM-Newton temperatures and calculate the projected mass, which is 
consistent with the XMM-Newton mass and a little lower than the weak lensing mass at rsoo- Thus, it seems that 
the temperature derived from XMM-Newton may be more reliable. 
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1. Introduction 

Galaxy clusters, the largest objects in the universe, origi- 
nate in the primordial density perturbations from cosmic 
gravitational collapse. They are used in a variety of ways 
to perform both cosmological and astrophysical studies. 
Modern astronomy satellites such as XMM-Newton and 
Chandra have high sensitivity and spatial resolution, and 
unprecedented results have been gained through detailed 
analysis of their data. However, there are discrepancies 
between the properties of galaxy clusters derived from 
Chandra and XMM-Newton data, such as gas tempera- 
ture and total mass. Nevalainen et al. (2010) examined 
the cross-calibration of the energy dependence and nor- 
malization of the effective area of Chandra and XMM- 
Newton, finding that the discrepancies of the 0.5 — 7.0 
keV band temperature measurements of galaxy clusters 
with EPIC '/XMM-Newton and ACIS/ Chandra could reach 
i~ 10 — 15% on average. 

Abell 1835 is a classical bright cluster with a big cool 
core (Allen et al. 1996). Its X-ray morphology (Schmidt 
et al. 2001; Smith et al. 2005) shows that it is an undis- 
turbed and relaxed cluster. It is also an optimal candidate 
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for a triaxial joint analysis via X-ray, SZ, and lensing tech- 
niques (Morandi et al. 2011). The temperature difference 
between the Chandra and XMM-Newton analysis also ex- 
ists in Abell 1835: ~ 12 keV from Chandra data (Schmidt 
et al. 2001) and ~ 7.6 keV from XMM-Newton data (Jia 
et al. 2004). After the point-spread function (PSF) cor- 
rection of XMM-Newton data, Wang et al. (2010) derived 
a temperature profile similar to that of Jia et al. (2004). 
Therefore, the difference is not due to the PSF effect. 

In addition, the masses of Abell 1835 derived from 
XMM-Newton and Chandra are different (Jia et al. 2004; 
Schmidt et al. 2002). Fortunately, the masses have been 
measured with the strong lensing method (Richard et al. 
2010) as well as with the weak lensing method (Zhang et 
al. 2008). Since gravitational lensing directly probes the 
cluster total mass without any strong assumptions about 
the equilibrium state of the cluster, the lensing mass is 
generally more reliable. 

Recently, the comparison of X-ray and gravitational 
lensing masses has been studied in detail by both ob- 
servational and simulated analyses. Generally, the X-ray 
mass is consistent with or lower than the gravitational 
lensing mass. Richard et al. (2010) showed that the ratio 
of strong lensing mass and X-ray projected mass within 
r < 250 kpc, Msl/Mx, was 1.3. It was also found that the 
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mass derived from the X-ray measurement is about half 
of the strong lensing mass in some clusters (e.g., A1689, 
Andersson & Madejski 2004, Lemze et al. 2008; PKS 0745- 
191, Chen et al. 2003). Moreover, Zhang et al. (2008) 
showed that the average ratio of the weak lensing mass to 
X-ray mass was 1.09 ± 0.08, while Mahdavi et al. (2008) 
demonstrated that M x /M W l is 1.03±0.07 and 0.78±0.09 
at r 2 5oo and r 50 o, respectively. N-body/hydrodynamical 
simulation work also estimated the ratio between X-ray 
and lensing masses, Mx/Mwl, which was 0.88 ±0.02 and 
0.75 ± 0.02 at r 500 in Meneghetti et al. (2010) and Rasia 
et al. (2012), respectively. 

In this work, we reanalyse the Chandra data of Abell 
1835 with the same deprojection technique as Jia et 
al. (2004) to ascertain if the differences in temperatures 
derived from Chandra and XMM-Newton data can be 
corrected by data analysis. Furthermore, we investigate 
the reasons for the temperature differences and establish 
which result is more reasonable. 

This paper is organized as follows: The Chandra obser- 
vations and data preparation are described in Section 2. 
Section 3 shows the basic spectra analysis, while Section 
4 presents the deprojected electron density profile and the 
mass profile. In Section 5, we discuss the reasons for the 
temperature differences derived from Chandra and XMM- 
Newton data. We present our conclusions in Section 6. 

Throughout this paper, we assume Hq — 70 km 
s _1 Mpc _1 , f^A = 0.7, il m = 0.3. One arcminute corre- 
sponds to 236.2 kpc at Abell 1835 redshift of 0.2523. The 
selected energy band is 0.5 keV- 7.0 keV. 

2. Observation and spectra extraction 

Abell 1835 was observed with Chandra on 25 August 2006 
for a total of 119.48 ksec (observation ID is 6880). The 
observation instrument was ACIS-I, and the observation 
model was VFAINT. We processed the Chandra data with 
CIAO 4.2 and CALDB 4.2.0. 

Background images were extracted from the stan- 
dard set of CTI-corrected ACIS blank sky images in the 
Chandra CALDB (Markevitch et al. 2003). To remove the 
particle background, we estimated the count rates in the 
hard energy band (10-12 keV) of Abell 1835 and the blank 
sky, and renormalized the blank sky. We used the tool 
LC_CLEAN in CIAO to scan the light curve of data for 
flares. The Good Time Interval (GTI) was about 117 ksec. 
Because Abell 1835 appears to be a relaxed cluster of 
galaxies, we assumed that the temperature structure of 
this cluster is spherically symmetric and applied a depro- 
jection technique. The deprojected spectrum of each shell 
is calculated by subtracting the contributions from all the 
outer shells (e.g., Nulsen & Boringer 1995; Matsushita et 
al. 2002). The detailed calculation procedures were de- 
scribed in Chen et al. (2003) and Jia et al. (2004, 2006). 

We divided the image of the cluster into seven annu- 
lar regions centered on the emission peak for the extrac- 
tion of spectra and used the outmost ring (8.33' — 10.42') 
to determine the local cosmic X-ray background (CXB). 



For each annular region, ancillary response files (ARFs) 
and response matrix files (RMFs) are generated using 
CIAO. The complete process is: a) select GTI using the 
light curve of data; b) subtract the point sources; c) sub- 
tract the background using blank sky data; d) extract the 
project spectra; e) extract the deprojected spectra. 

For the purpose of comparing, we also reprocessed 
the XMM-Newton data of Abell 1835 (observation ID: 
0551830201) in ACDM cosmology as done in Jia et al. 
(2004) using SAS 11.0.0. 

3. Spectral analysis 

For the spectral analysis, we used the plasma emission 
model MEKAL (Mewe et al. 1985, 1986; Kaastra 1992; 
Liedahl et al. 1995) and WABS model (Morrisson & 
McCammon 1983). To fit the spectra, XSPEC version 
12.6.0 (Arnaud 1996) is used, and the model is 

Modeh = Wahs(n H ) x Mekal(T, z, A,norm). (1) 

We fixed the redshift z to 0.2523 and n# to the Galactic 
absorption 2.24 xlO 20 cm~ 2 (Dickey & Lockman 1990). 
The fitting results are listed in Table 1. 

The abundance is higher in the cluster center, which 
is understandable because the excess metal in the cluster 
center is produced in the cD galaxy and ejected into the 
ICM (Makishima et al. 2001, Xiang et al. 2009). From 
the deprojected temperature profile (squares in Fig.l), we 
see that the temperature decreases towards the center, 
which may be ascribed to the gas cooling. Here we find 
the temperature profile is fitted well with the formula 

T(r)^ae-^- b ^/ 2c \ (2) 

The best-fit parameters arc a = 11.29 keV, b = 2.78 ar- 
cmin, and c = 2.14 arcmin. The best-fit profile is shown 
as a solid line in Fig.l. 

To establish a comparison, we also obtained the tem- 
peratures from the projected spectra of Chandra, plotted 
as diamonds in Fig.l. The average temperature excluding 
core (r < 0.75') is ^9.6 keV, which does not differ much 
from that of Markevitch (2002). After deprojection, the 
temperatures do not differ much from the projected tem- 
peratures. But they are still much higher than our new 
result of XMM-Newton, ~7.33 keV, which is consistent 
with that of Majerowicz et al. (2002) and Jia et al. (2004). 
Therefore, even when the same method of data analysis 
is used, the temperature differences derived from Chandra 
and XMM-Newton data still exist. 

4. Mass analysis 

In our mass model, we measure the spatially resolved ra- 
dial temperature distribution from the deprojected spec- 
tra. A double-/? model is adopted to fit the ICM density 
distribution. 
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Table 1. Best-fit free parameters of Abell 1835: temperature T, abundance A, and normalized constant norm. The errors 
represent a confidence level of 90%. 
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Fig. 1. Temperature profiles of Abell 1835 with a confidence 
level of 90%. Squares: deprojected temperatures of Chandra; 
diamonds: projected temperatures of Chandra; triangles: de- 
projected temperatures of XMM-Newton. The solid line is the 
best-fit profile of Chandra deprojected temperatures. 



4.1. Electron density 

For calculating the electron density profile, we divided the 
cluster into 19 annular regions. We calculated the depro- 
jected photon counts in each shell. By using the depro- 
jected abundance and the deprojected temperature pro- 
file, we could estimate the normalized constant of each 
region, norm. Then we derive the deprojected electron 
density n e of each region from Eq(3), shown as squares in 
Fig.2, 



10- 



norm — 



4ir[D A (l + z)} 2 



n e n H dV, 



(3) 



where Da is the angular size distance to the source 
(cm) and n ei nn (cm -3 ) are the electron and hydrogen 



densities, respectively. We fitted the electron density pro- 
file with the double-/? model (Chen et al. 2003) 



n e (r) = n i 



1+ — 



r 
r c i 



+n 2 



•(4) 



The best-fit parameters of Chandra arc noi = 0.041 ± 
0.002 cm -3 , /3i = 0.597 ± 0.006, r cl = 0.325 ± 
0.013 arcmin, n a2 = 0.114 ± 0.003 cm~ 3 , /3 2 = 12.81 ± 
3.080, r c2 = 0.730 ± 0.094 arcmin, xl ed = 23.48, d.o.f. = 
13. 

The best-fit profiles of n e from Chandra and XMM- 
Newton data are depicted in Fig.2, which shows that the 
electron density derived from Chandra is higher than that 
of XMM-Newton. 



4.2. Mass calculation 

With the assumption of spherical symmetry and hydro- 
static equilibrium, the total mass of cluster within radius r 
can be determined when the radial profiles of the gas den- 
sity and temperature are known. We calculated the gravi- 
tational mass of Abell 1835 with the hydrostatic equation 
(Fabricant et al. 1980) 



M tot {< r) 



k B Tr 2 r d(lnn e ) d(lnT) 1 



Gfim p dr 



dr 



(5) 



where Ub is the Boltzmann constant, G is the gravitational 
constant, \i is the mean molecular weight of the gas in 
units of m p , and m p is the proton mass. For a fully ionized 
gas with a standard cosmic abundance, a suitable value is 
fj, = 0.60. 

Using the best-fit profiles of the electron density n e (r) 
and the deprojected temperature T(r), we could obtain 
the total mass profile, as shown in Fig. 3. The total masses 
within the radius of 6' are M to t-Chandra — 1-36 ± 0.42 x 
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Fig. 2. Top panel: electron density profiles of Abell 1835 from 
Chandra (squares) and XMM-Newton (triangles). The error 
bars are at a 68% confidence level. The solid line is the best- 
fit profile of Chandra results with a double-/? model, and the 
dashed line for XMM-Newton results. Bottom panel: the ratio 
between the electron density profiles of Chandra and XMM- 
Newton. 



M SL = 2.83±0.41 x 10 14 M Q , based on seven multiply im- 
aged systems under ACDM cosmology. We estimated the 
X-ray projected mass inside the radius M arc -chandra = 
1.85 ± 0.20 x 10 14 M Q , M arc _xMM = 1-40 ± 0.15 x 10 14 
M©. The results, shown in Fig. 4, are lower than the strong 
lcnsing mass derived by Richard et al. (2010). The differ- 
ence may be due to the assumption of spherical symmetry 
and hydrostatic equilibrium in the X-ray mass calculation 
(Gavazzi 2005). 

The mass analysis of strong lensing cluster may be 
affected by orientation biases (Meneghetti et al . 2010b, 
2011; Hennawi et al. 2007; Oguri & Blandford 2009; Zitrin 
et al. 2011). For Abell 1835, Corlcss et al. (2009) showed 
that the assumption of spherical symmetry is less prob- 
lematic and that the orientation bias is weaker for this 
halo because of its smaller Einstein radius. Thus, the mass 
estimates of Abell 1835 mentioned above are reasonable. 
The radio plasma may provide some additional pressure 
to support the X-ray gas, so the assumption of hydrostatic 
equilibrium may be inaccurate. 



10 15 M© and M tot _ XMM = 0.84 ± 0.10 x 10 15 M . This 
shows that the Chandra result is much higher than the 
XMM-Newton result. At the same time, the virial mass of 
Chandra is M 200 = 2.0±g;| x 10 15 M© with r 200 = 1-99 
Mpc, which is consistent with the result of Schmidt et al. 
(2001) when using the same cosmology model. 



14 — 



12 — 




Radius (arcmin) 

Fig. 3. Total mass profiles of Abell 1835 from XMM-Newton 
(triangles), Chandra (squares) and a method in subsection 5.2 
(circles) fitting Chanda spectra with XMM-Newton tempera- 
tures. 



4.3. Mass within the optical lensing arc 

Allen et al. (1996) announced that there is a lensing 
arc inside the cluster center. Richard et al. (2010) accu- 
rately recalculated the mass inside the radius of 250 kpc, 
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Fig. 4. Projected mass profiles within 250 kpc of Abell 1835 
from Chandra (squares) and XMM-Newton (triangles). The 
dashed line presents the NFW profile of strong lensing mass 
of Richard et al. (2010). 



4.4. Gas mass fraction 

In galaxy clusters, gas , which has a temperature of a 
few keV, is an important component. From the electron 
density, we can calculate gas mass and then the gas mass 
fraction, defined as f gas {r) = M gas (r)/M tota i{r). 

Fig. 5 shows the gas mass fraction of Abell 1835, which 
drops clearly in the region r < 0.6', while keeping approx- 
imately constant in the outer region. This indicates that 
the dark matter is more condensed than the gas in the 
central region of the cluster. 
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Fig. 5. Gas mass fraction profiles of Abell 1835 from Chandra 
(squares) and XMM-Newton (triangles). 

5. Discussion 

Although we used the same method to analyse the data 
of Abell 1835 observed by XMM-Newton and Chandra, 
the temperature discrepancy still exists. The average tem- 
peratures of Abell 1835 (excluding cool core) derived by 
Chandra and XMM-Newton are 9.6 keV and 7.33 keV, re- 
spectively, in the following, we discuss the possible reasons 
for this. 

5.1. PS F scatter 

The effect of PSF of pn/XMM-Newton and MOS/XMM- 
Newton is larger than that of ACIS/ Chandra. PSF may 
scatter a significant fraction of the adjacent regions' emis- 
sion into the regions that we analyze. 

Applying a combined direct demodulation and depro- 
jection technique to subtract the effects of PSF, Wang 
et al. (2010) found that the central electron density in- 
creases by 30%, while the temperature profile is similar to 
the result using only a deprojection method. They also es- 
timated that the effects of PSF were not important for the 
temperature profile. Nevalainen et al. (2010) pointed out 
that PSF scattering may cause differences in temperature 
measurements with different instruments of a maximum 
of 2%, which is negligible compared to the statistical un- 
certainties with the minimum width of rings set to 1.5'. 
Therefore, PSF is not a primary reason for the tempera- 
ture discrepancy. 

5.2. Discrepancy between the temperatures 

Since the effects of PSF and deprojection technology on 
temperature are not the primary reasons for the tempera- 
ture discrepancy, the difference of the calibration of XMM- 
Newton and Chandra may lead to the temperature discrep- 
ancy. Nevalainen et al. (2010) found that the difference 
of the calibration of these two instruments is about 10- 
15% on average and the maximum difference is ~ 25%. 



Therefore, it may be the main factor of the temperature 
discrepancy of Abell 1835. 

In order to contrast the temperatures of Abell 1835 
determined by Chandra and XMM-Newton, we tried to 
fit Chandra and XMM-Newton spectra with temperatures 
fixed to the other one's temperatures. Because the tem- 
peratures of XMM-Newton and Chandra are significantly 
different in annulus 0.75' — 3.33', we combined and re- 
binned the three spectra of 0.75' — 3.33'. The fitting re- 
sults are listed in Table 2. It seems that Chandra spectra 
can be fitted well with XMM-Newton temperature, while 
XMM-Newton spectra can not fit well with Chandra tem- 
perature. As a result, we think that XMM-Newton has a 
stronger restriction on temperature. 

Furthermore, we calculated the total mass with 
XMM-Newton temperatures and Chandra electron den- 
sity based on XMM-Newton temperatures (hereafter 
Mchandra(T=T XMM )), which is 0.74±0.05 x 1O 15 M within 
the radius of 6.0' (circles in Fig. 3). It is similar to the total 
mass derived from XMM-Newton data. 

To compare with the weak lensing mass, we also cal- 
culated the projected masses at r25oo, ^1000 an d rsoo- The 
results are listed in Table 3, which shows that the weak 
lensing mass is lower than Mchandra, and a little higher 
than Mxmm and M chandra ( T=TxMM) at r 500 . 

Because of the non-hydrostatic state and non- 
equilibrium processes in clusters, the X-ray masses may 
be underestimated (e.g., Rasia et al. 2004; Piffaretti et al. 
2004; Nagai et al. 2007). In Table 4, we list the compar- 
ison of the X-ray mass and the weak lensing mass from 
both the observational and simulated analyses. It shows 
that the X-ray mass is always consistent with, or a little 
smaller than, the weak lensing mass (Zhang et al. 2008; 
Mahdavi et al. 2008; Zhang et al. 2010; Meneghetti et al. 
2010; Rasia et al. 2012). Consequently, our Mxmm and 



M, 



Chandra(T=T X MM) 



are more suitable. 



With XMM-Newton temperatures, we can get more 
reliable masses from XMM-Newton and Chandra data. 
Because XMM-Newton has a stronger restriction on tem- 
perature, the temperatures derived from XMM-Newton 
data may be more reliable. 

6. Conclusion 

We have presented a detailed analysis of about 117 ksec 
of Chandra observations on the galaxy cluster Abell 1835. 
Through the deprojected spectra analysis, we derived the 
deprojected temperatures, which do not differ much from 
the temperatures without deprojection and are still much 
higher than those of XMM-Newton. 

The total mass within the radius of 6', M to t-Chandra — 
1.36±0.42xl0 15 M Q , is much higher than the total mass 
of Mtot-XM m = 0.84±0.10xl0 15 M© derived from XMM- 
Newton data. The difference of total mass is due to the 
discrepancy of the temperatures derived from Chandra 
and XMM-Newton. We also calculated the projected mass 
within the optical lensing arc, which is lower than the 
strong lensing mass derived by Richard et al. (2010). 
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Table 2. Best-fit parameters of 0.75' — 3.33': temperature, abundance, and normalized constant norm. The errors represent a 
confidence level of 90%. 



Parameter 


Chandra 


Chandra 
(T = Txmm) 


XMM-Newton 


XMM-Newton 

(T — Tchandra) 


Temperature (keV) 


10.44+°-" 

— 0. 1 2, 


7.73 (fix) 


7.73+°S 

— U.o4 


10.44 (fix) 


Abundance (solar) 


n 97 +0.09 
u '^' -0.09 


0.27 (fix) 


n S1+ 09 

U-Ol_g 0g 


0.31 (fix) 


norm (10~ 3 cm~ 5 ) 


4 1 7 +0.07 
' -0.07 


q - lo -0 .08 


4 Q4+0.O6 


4 Q4+0.O6 


Xred/d-O.f. 


0.86/191 


1.24/193 


0.98/165 


1.51/167 


Probability 


0.91 


0.013 


0.56 


2.24 x 10~ 5 



Table 3. Projected masses of Abell 1835. The weak lensing masses are reported in Zhang et al. (2010). r2soo, »*iooo and rsoo 
are determined from the weak lensing analysis. 



R 


Mchandra 

10 14 M 


MxMM — Newton 
10 14 M 


Mchandra(T=T KMM ) 

10 14 M 


M w i 
10 14 M 


5*2500 


4.42 ± 0.52 


3.15 ±0.62 


3.28 ± 0.40 


2.88 ±0.58 


**1000 


7.11 ± 1.86 


5.80 ± 1.42 


5.28 ± 1.39 


6.15 ±0.95 


7*500 


13.80 ±3.91 


8.99 ± 3.60 


8.19 ±3.21 


9.65 ±1.70 



Table 4. Ratio between X-ray and weak lensing masses. 

Sample M x /M W l 
-R2500 -fi+ooo R500 

Observation Zhang et al. 2010 (unrelaxed) 1.00 ± 0.07 0.97 ± 0.05 0.99 ± 0.07 

Zhang et al. 2010 (relaxed) 1.04 ± 0.08 0.96 ± 0.05 0.91 ± 0.06 

Mahdavi et al. 2008 1.03 ± 0.07 0.90 ± 0.09 0.78± 0.09 

Simulation Rasia et al. 2012 0.83 ± 0.02 0.80 ± 0.02 0.75 ±0.02 

Meneghetti et al. 2010 0.90 ±0.04 0.86 ± 0.02 0.88±0.02 



After dcprojection and PSF correction (Wang et al. 
2010), the temperature and total mass of Abell 1835 
resulting from Chandra are still different from those of 
XMM-Newton. These differences may result from the cal- 
ibration of the two instruments. We find that XMM- 
Newton has a stronger restriction on temperature. And 
with XMM-Newton temperatures, the projected mass 
from Chandra data is lower than the weak lensing mass 
and consistent with other observational analyses. For these 



reasons, the temperatures obtained from XMM-Newton 
may be more reliable. 
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